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Abstract 

The authors measure the persistence and contour lengths of DNA-psoralen complexes, as a 
function of psoralen concentration, for intercalated and crosslinked complexes. In both cases, 
the persistence length monotonically increases until a certain critical concentration is reached, 
above which it abruptly decreases and remains approximately constant. The contour length of 
the complexes exhibits no such discontinuous behavior. By fitting the relative increase of the 
contour length to the neighbor exclusion model, we obtain the exclusion number and the intrinsic 
intercalating constant of the interaction. Ab initio calculations are employed in order to provide 
an atomistic picture of these experimental findings. 

PACS numbers: 87.80.Cc, 87.14.gk, 87.15.La 
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DNA interactions with ligands such as drugs or proteins have been extensively studied 



and reviewed by many authors in the past years 



-[8|. Some of these drugs, like daunomycin 
3|, |6j. Psoralen, by contrast, 



and ethidium bromide (EtBr), exhibit intercalative binding jl, 
presents other forms of linkage to the DNA bases. In fact, it is well known that this drug 
can absorb photons and form covalent bonds with the DNA bases, when the sample is 
illuminated with ultraviolet-A (UVA) light (A = 320-400 nm) 4|. Psoralens are compounds 
from the family of furocoumarins, broadly used in the medical treatment of various skin 
diseases like psoriasis, vitiligo, and some other kinds of dermatitis {9]. 




FIG. 1: Psol-DNA fragment with five base CG pairs and two intercalated psoralens obtained from 
our ab initio DFT calculations. 



Although psoralens have been extensively used in medicine for many years, the interac- 
tion of the drug with DNA is not yet completely understood. In this work, we study the 
intercalative binding and crosslink formation in DNA-psoralen complexes (Psol-DNA), in or- 
der to clarify the quantitative aspects of the interaction between psoralen and DNA. We are 
interested in probing the modifications of the mechanical properties of DNA molecules when 
interacting with psoralen. The elasticity and size changes of these molecules, as a function of 
the total psoralen concentration (Cp), are measured by performing single- molecule stretch- 
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ing experiments, which give information about persistence (£) and contour (L) lengths of 
Psol-DNA. The total psoralen concentration is the sum of the intercalated concentration (C,) 
and the free concentration in solution of psoralen. We also define the intercalated psoralen 
fraction r=Ci/Cb p , where Cb p is the DNA base-pair concentration, fixed in all experiments 
at 11 fiM. In order to elucidate the microscopic mechanisms behind the psolaren-DNA in- 
teraction, we also perform first principles calculations for Psol-DNA fragments (an example 
is shown in Fig. [Q to determine the behavior of the Psol-DNA stiffness (which is related 
to the Young modulus) as a function of Cp, and also to obtain the maximum fraction of 
intercalated psoralen into DNA, i.e., the limiting value of Cp at which a Psol-DNA complex 
is still structurally stable. 

Experimentally, we observe that as Cp increases the Psol-DNA persistence length initially 
increases and then undergoes an abrupt transition to a lower value at a critical total psoralen 
concentration (Cp lt ), a behavior we also observe in UV-light illuminated complexes. In 
the latter case, the Psol-DNA complex attains higher persistence length values for Cp < 
Cp rit , which is indicative of crosslinking, where the psoralen binds covalently preferably to 
thymine bases in opposite DNA strands, upon illumination. As discussed below, the ab initio 
calculations are in qualitative agreement with these experimental findings, enabling us to 
analyze the changes in the molecular structure that are responsible for the experimentally 
observed behavior. 

The intercalated complexes are obtained by simply waiting for the psoralen to interca- 
late the A-DNA in the sample. The crosslinked complexes are obtained by illuminating the 
sample with a UVA mercury lamp. We use an optical tweezers to trap a polystyrene bead 
attached to the DNA molecule, while pulling the microscope stage with a controlled veloc- 
ity. We then make force x extension curves, and use the approximate expression derived by 
Marko and Siggia [h]] to obtain the persistence and contour lengths of the bare DNA and 
Psol-DNA complexes. The persistence length £ = YI/kT (Y is the Young modulus, I is the 
geometrical moment of inertia of the polymer, and kT is the thermal energy) is a mesoscopic 
quantity which depends on the local elasticity of DNA. Its value is £ ~ 50 nm for bare DNA 
under physiological conditions. In our experiments, the forces applied to stretch the DNA 
are within the entropic force regime (forces < 3 pN), thus avoiding externally distorting the 
structure and shifting the DNA/psoralen chemical equilibrium (enthalpic effects). The de- 
tails about our experimental setup, experimental procedure, and optical-tweezers calibration 
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FIG. 2: Persistence length £ (in nm) as a function of total psoralen concentration Cp (in /iM) for 
Psol-DNA complexes. Circles: intercalated Psol-DNA; Squares: crosslinked Psol-DNA. The DNA 
base-pair concentration used was Cb p =ll fiM. 



can be found in Refs. 



file 



In an attempt to understand the behavior of £ as a function of Cp, we perform first prin- 
ciples calculations of the structural stability and stiffness changes produced by intercalation 
of psoralen into fragments of dry poly(dG)-poly(dC) in its acid form. We investigate DNA 
models with fragments consisting of four and five guanine-citosine (CG) pairs. A five-base- 
pair fragment with two intercalated psoralens is shown in Fi g. [1 Our first principles density 
functional theory (DFT) simulations use the SIESTA code [15[. Exchange and correlation 
effects are described within the generalized gradient approximation [l^ . A double-^ basis set 
of atomic orbitals of finite range is used, with polarization functions added for phosphorus 
and for the atoms involved in the hydrogen brid ges. This methodology has been shown to 



provide a good description of the DNA structure 



17 



191 ] . We treat isolated finite fragments 



in a periodic supercell with large vacuum regions. For the DNA bases at the edges of the 
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fragments, one atom was held fixed in the same position as in the infinite periodic DNA [17J. 
Broken bonds at the edges are saturated with hydrogen atoms. Full geometry relaxations 
were performed (forces < 0.04 eV/A). The stability of Psol-DNA complexes is addressed by 
computing their formation energies defined as the difference between the total energy of the 
Psol-DNA complex and the sum of the total energies of the bare-DNA fragment and the 
isolated psoralen molecule, after full relaxation of the initial geometries. 
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FIG. 3: Total psoralen concentration Cp as a function of relative increase of the contour length 
for Psol-DNA intercalated complexes. for psoralen molecules equals the intercalated psoralen 
fraction r. Dots are the experimental data and the dashed line is a fit using Eq. 1. 



In Fig. |2l we plot the measured values of £ as a function of Cp, for both intercalated 
and crosslinked Psol-DNA complexes. The figure shows that £ increases monotonically 



until the critical concentration C 



crit 



13.4 /iM is reached. For Cp > Cp , £ decays 



abruptly to around 50 nm, and remains approximately constant at this value, at least for 
the concentration range in our experiments. Note that Cp lt is the same for both intercalated 
and crosslinked complexes. A similar abrupt transition of the persistence length was recently 
reported by two of us for other intercalating drugs (daunomycin and ethidium bromide) 6]. 
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The fact that Psol-DNA complexes exhibit this same behavior is a strong evidence that this 
transition is of general character for intercalating molecules, at least in the low-force regime 
(< 3 pN) used in our experiments. This abrupt transition may be caused by local formation 



of denaturing bubbles, as proposed in Refs. 



121 ] and 13] , which softens the structure of the 



Psol-DNA complex. In Ref. 12[ large scale molecular dynamics calculations indicate the 
formation of denaturing bubbles when pure DNA is stretched above 65 pN forces. Since we 
use very small stretching forces it is likely that the abrupt transition we observe is caused 
by some intrinsic DNA structural changes induced by the intercalating drug. 

In Fig. O we show the relative increase of the contour length, = (L — L Q )/L , as a 
function of Cp for the same intercalated complexes shown in Fig. [2J where L Q is the contour 
length for Cp = 0. Note that the contour length does not exhibit any abrupt changes. 
The dashed line in this figure is a fit using Eq.l below, derived from the neighbor exclusion 
model (NEM) Q (see Ref. |s| for details). Considering that each intercalated psoralen 
molecule increases the length of DNA by the base-pair distance 0.34 nm {4], is equal to 
the intercalated psoralen fraction r. In this case the NEM gives 

eq-ne + er 1 

° P ~ Cbp ° + K t (l-n&r ' { ) 

where n is the exclusion number and Ki is the intrinsic intercalating constant. From the 
fitting we determine the parameters n = 1.43 ± 0.13 and K\ = (8.8 ± 2.4) x 10 4 M _1 , 
for intercalative binding of psoralen. The large error bars for are due to the fact that 
different DNAs were used, and there is a natural length distribution for A-DNA. The abrupt 
transition for £ occurs around =0.38, while the maximum psoralen intercalated fraction 
(1/n) is between 0.64 and 0.77. 

Figure H] shows the behavior of the stiffness, which is proportional to the Young modulus 
(for fixed DNA length), as a function of r, computed using the DFT methodology. Note 
that the stiffness increases as r increases, and an abrupt decrease is observed for r between 
0.40 and 0.50, in reasonably good agreement with the data for the persistence length of 
Fig |2j where this transition occurs around = r =0.38. Our calculations indicate that 
before the transition two psoralens are well intercalated into the DNA structure, while after 
the transition one of the psoralens binds with the citosine, and the hydrogen bond between 
the two bases is broken, indicating that local denaturation may indeed be the mechanism 
behind the abrupt change in stiffness and consequently in the persistence length. Moreover, 
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FIG. 4: Stiffness as a function of intercalated psoralen fraction (r), from ab initio DFT calculations. 
Results for r = 0, 0.25, and 0.50 were obtained using a four-base-pair model. Results for r = 
0, 0.20, 0.40, and 0.60 were obtained using a five-base-pair model. 

our calculations also indicate that when r increases above 0.60, the intercalated structure 
is unstable, which sets an upper bound on the maximum fraction of intercalated psoralen 
within this theoretical model, in good agreement with the lower-limit value of 0.64 measured 
experimentally. Our interpretation for the microscopic nature of the transition rests on the 
good qualitative and quantitative agreement between theory and experiment. We must, 
however, not refrain from stating the limitations of our model system, since calculations 
in such small fragments are strongly affected by boundary effects. This is reflected in the 
large values we obtain for the theoretical stiffness, which are also due to the fact that such 
small fragments are much stiffer than the much longer DNA chains in the experiments (The 
stiffness x length we compute for bare DNA using the four-base-pair model is about 20% 
higher than that of the five-base-pair model.) 

In conclusion, the authors observe experimentally that the persistence length of DNA- 
psoralen complexes increases with psoralen concentration up to a critical concentration, 
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decreasing abruptly and remaining approximately constant above this concentration, a be- 
havior also observed in UV-light illuminated crosslinked complexes. Ab initio calculations 
show a similar behavior for the stiffness of psoralen-DNA complexes, and also indicate that 
local denaturation and binding between psoralen and a DNA base may be the mechanism 
behind this persistence length transition. 

The authors acknowledge support from the European Grant EU (FP6-029192), and from 
the Brazilian agencies: FAPEMIG, CAPES, CNPq, Instituto do Milenio de Nanotecnolo- 
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